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SU Aurigae is a classical T Tauri star of spectral type G2. This star was 
one of the scientific targets of the MUSICOS 96 multi-site campaign that 
provided a wealth of high resolution cross-dispersed spectral data, with a 
good continuous time coverage (Unruh et al. 1998a). We present the re- 
sults of the analysis of the complex circumstellar environment of this star, 
with particular regard to magnet ospheric models, in which the accretion 
(S| ■ from the disk is channelled onto the star along magnetic field lines. The 

signatures of modulated outflows and mass accretion events are present 
■ in the spectra, as well as transient spectral features. We computed auto- 

correlation and cross-correlation functions to better investigate the source 
4^ ■ of the profiles' variability. The comparison of the profiles of different spec- 

^! tral lines allows us to study the footprints of events effectively observed 

O ' a l different distances from the stellar surface. 

• • 

Introduction 

X 

SU Aurigae is a G2 classical T Tauri star (cTTS). It has a projected rotational 
velocity of ~ 66 km s _1 and a rotational period of about 3 days. So far the best 
period estimates come from the analysis of the spectral variations of the Balmer 
lines. Giampapa et al. (1993) found a period of 2.98 ± 0.4 days in the blue wing 
of the Ha spectral line, later confirmed also in the red wing of H/3 by Johns & 
Basri (1995b). Petrov et al. (1996) reported a period of 3.031 ± 0.003 days mea- 
sured on the redshifted variations of the Balmer lines. Our data set suggests a 
somewhat shorter rotational period of ~ 2.8 days (Unruh et al. 1998a,b). Using 
the Ha, H/3, Na I D and He I D3 profiles, we try to disentangle the two main 
components that are believed to be present, namely accretion and wind signa- 
tures. From equivalent width measurements of fitted components in Ha and H/3, 
Johns & Basri (1995b) advocated that the accretion (Hj3 red absorption feature) 
and wind (Ha blue absorption feature) signatures are out of phase in SU Aur, in 
what they call the misaligned "egg-beater" model or the oblique rotator, as it is 
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also known. This model is a generalization of the magneto-centrifugally driven 
flow model of Shu et al. (1994). According to this model, a stellar dipole field 
truncates the accretion disk close to the co-rotation radius. At the truncation 
point, the ionized disk material is loaded either onto inner closed magnetic field 
lines, accreting thereby onto the stellar photosphere, or it is loaded onto outer 
open magnetic field lines that can drive a disk- wind flow. In this context, we 
analyse H/3, Na I D and He I D3 data sets in terms of cross-correlation tech- 
niques (expanding on the work of Johns & Basri (1995a) for the Ho: spectral 
line). Also intriguing in our data set are two drifting features appearing in Ha 
and H/3 and, in one case, also in the Na I D lines. 



2. Description of the Data Set 

A more detailed description of our data set can be found in Unruh et al. (1998a), 
but we highlight here the most important characteristics. This data set was ob- 
tained in November 1996 during the MUSICOS 96 multi-site, multi-wavelength 
campaign, that involved five different observatories: Isaac Newton Telescope 
(INT, La Palma), Observatoire de Haute Provence (OHP, France), McDonald 
Observatory (USA), Xinglong Observatory (China) and Canada- France-Hawaii 
Telescope (CFHT, Hawaii). We obtained 126 echelle spectra spanning 8.5 days, 
with 3.5 days covered almost continuously. Ha, Na I D and He I D3 were ob- 
served at all five observatories while H/3 could only be observed from the INT, 
OHP and CFHT. Our data set differs considerably from previous SU Aur data 
sets, because those had typically one to two spectra a night, even though over 
longer time spans. Our data set, by contrast, has a finer time sampling, so we 
are less limited by the loss of coherency of the observed phenomena. We can 
hence study the timescales over which the variations of the different spectral 
lines are related. 

3. Cross-Correlation Analysis 

The broad spectral coverage of our data set allows the comparison of lines that, 
due to their different ionization potentials, probe different parts of the circum- 
stellar material. We compute the cross-correlation function of pairs of spectral 
lines as a function of the time lag (At). The data sets are interpolated to ac- 
count for the unequal spacing (White & Peterson 1994). In simple terms each 
velocity bin in one spectral line (at time t) is correlated with all the velocity 
bins of the other spectral lines (interpolated at time t+At) and this gives, for 
each pair of lines, an intensity map like the examples shown in Fig. |l|. These 
maps provide a wealth of information , even though not being necessarily the 
easiest to interpret. Starting with the auto-correlation function of H/3 (see top 
row in Fig. |l|): 

• The red wing of H/3 auto-correlates very well over the interval [50:200] km s~ 

• The red wing of H/3 [50:200] km s _1 is anti-correlated with the blue wing 
of H/? [-150:0] km s _i . 
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Figure 1. Cross-correlation intensity maps for H/3, H/3 versus Na I D 
and H/? versus He I D3. Each row represents intensity maps for one pair 
of lines, for 3 time lags, At = 0, At ~ P/2, At ~ P. Yellow represents 
high correlation coefficients and black strong anti-correlations. The 
intensity (color) scale is not the same in all maps. In each map, the 
time lag (At) is indicated as well as the maximum/minimum correlation 
coefficients found. 
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• The red and blue wings are correlated albeit over a smaller velocity range 
for a time lag of ~ 1.5 days, approximately half of the detected period. 
For the same time lag the red wing is anti-correlated with itself. 

• With a time lag of ~ 2.8 days, the period we detected, we recover almost 
the initial cross-correlation map, due to periodic variations on both wings. 

Thus the blue and red wings of H/3 seem indeed to be 180° out of phase. For 
the cross-correlation of H/3 and Na I D a similar analysis was performed. The 
resulting maps are shown in the middle row in Fig. p]. The map at At=0 shows 
the same features as the H/3 cross-correlation map, although over a narrower 
velocity range. This suggests that the two Na I resonance lines behave very 
similar to the H/3 line. In particular we find: 

• the red wings of the two Na I D lines [25:150] km s _1 correlate well with 
the red wing of H/3 [50:200] km s" 1 . 

• With a time lag of about half of the detected period, these parts of the 
profiles become anti-correlated. 

• There is a weak correlation between the blue wings of Na I D [-100:0] km s _1 
and H/3 [-150:0] km s" 1 . 

• The red wings of Na I D [25:150] km s _1 are anti-correlated with the blue 
wing of H/3 [-150:0] km s" 1 . 

Therefore the red wings of Na I D and H/3 vary approximately in phase. The 
case of the cross-correlation between H/3 and He I D3 is slightly different, as the 
whole He I D3 seems to vary in a concerted way (see bottom row in Fig. D. We 
find that: 

• The whole He I D3 profile [-100:150] km s _1 correlates well with the red 
wins of H/3 [25:220] km s _1 ; in these regions the 2.8 day periodicity was 
detected. 

• With a time lag of about half of the detected period the correlation diagram 
shows these lines as anti-correlated. 

In conclusion, we find that the blue and red wings of H/? are approximately 180° 
out of phase, while the red wings of Na I D and H/3 vary in phase. Finally, the 
whole profile of He I D3 varies in phase with the red wing of H/3. The time 
lags here indicated have to be taken with caution, as no error estimates were 
performed. 

4. Oblique "Egg-beater" Model Predictions 

Johns & Basri (1995b) found that the blue wing of Ha and the red wing of 
H/3 were approximately 180° out of phase. They proposed that this behaviour 
could be explained by an oblique "egg-beater" model (see Fig. 15 of Johns & 
Basri 1995b) that predicts that the signatures of mass accretion and disk winds 
should be rotationally modulated and 180° out of phase, since at different phases 
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the visible wind flow or funnel flow are favoured. From the H/3 auto-correlation 
analysis we indeed have confirmation that the red wing and the blue wing ab- 
sorptions vary in this way, meaning that the accretion and disk- wind signatures 
are in phase opposition. Also favouring this model is the fact that the variance 
profile of H/3 shows in the red wing a preferred projected velocity (as expected 
from a funnel geometry); this is not seen in the blue wing, as this part of the 
profile originates in a less "collimated" disk- wind. The red wings of the Na I D 
lines vary in phase with the red wing of H/3, a clear sign that they are also 
formed in accretion flows. Futhermore, the red wing of H/3 varies in phase with 
the whole profile of He I D3. This can be taken as an indication that the helium 
line excitation is also related to the accretion but not to the wind in SU Aur. 
He I D3 excitation needs very high temperatures (or densities), therefore it can 
not be simply formed in the accreting flows. If the source of the excitation is 
the energy released at the "foot-points" of the accretion streams, then one could 
expect that the H/3 red wing (accretion) and the He I D3 line vary approximately 
in phase, as we have found. 



5. Transient Features 

A drifting feature appears between the eighth and eleventh days of our observa- 
tions in the profiles of Ho;, H/3 and Na I D (see Fig. 1 of Unruh et al. 1998b). By 
analysing the Na I D profiles, it is clear that this is an absorption component. 
As this part of the profile correlates well with the corresponding parts in the 
Balmer profiles, we treat it as the same drifting absorption feature, present in 
all three spectral lines. It is important to point out that no clear periodicity 
was detected in this part of the profiles for any of the lines. We have measured 
the velocity position of the minima of the line profiles in the relevant velocity 
interval and in Fig. [2] we show their velocity evolution. It seems to indicate that 
v(Hck) < v(H/3) < v(Na I D). However, the velocity measurements for H/3 have 
to be treated with caution as the feature can not be identified unambiguously 
in this line. Instead, we observe that the blueshifted absorption gets broader 
and deeper. The velocity difference (between the lines) probably traces material 
accelerating outwards. It should be stressed that approximately 1 day after the 
decay of the feature described above, another transient feature was detected in 
Ha and H/3 but no counterpart was seen in the Na I D lines. 

6. Conclusions 

The present analysis supports the so-called oblique or misaligned "egg-beater" 
model (Johns & Basri 1995b). Our data set covers 8.5 days with several spectra 
each night, and at least 3.5 days of almost continuous coverage, allowing us to 
compare different (parts of the) line profiles and study over which time scales are 
they related. The temporal relations we have found between the spectral lines 
help to clarify which parts of the line profiles variations originate from which 
component of the circumstellar environment of SU Aur: the red wing absorption 
components of the H/3 and Na I D lines form in the accretion funnels, He I D3 
forms in the highly energetic region at the "foot-points" of these columns, and 
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Figure 2. Velocity evolution of the drifting absorption feature seen 
in the Ha, H/3 and Na I Dl line profiles. For Na I D2, the velocity 
position of the feature is similar to Na I Dl and we do not represent it 
here for clarity's sake. 



the slightly blueshifted absorption feature in the blue wing of H/3 forms in a 
disk-wind flow. Also present in the data set are transient features not fully 
understood, possibly the signature of an outwardly accelerating (stellar) wind 
component. 
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